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ABSTRACT 

Context. Stellar limb darkening, I(p. = cos 6), is an important constraint for microlensing, eclipsing binary, planetary transit, and 
interferometric observations, but is generally treated as a parameterized curve, such as a linear-plus-square-root law. Many analyses 
assume limb-darkening coefficients computed from model stellar atmospheres. However, previous studies, using from plane- 
parallel models, have found that fits to the flux-normalized curves pass through a fixed point, a common fi location on the stellar disk, 
for all values of 7^ , log g and wavelength. 

Aims. We study this fixed /i-point to determine if it is a property of the model stellar atmospheres or a property of the limb-darkening 
laws. Furthermore, we use this limb-darkening law as a tool to probe properties of stellar atmospheres for comparison to limb- 
darkening observations. 

Methods. Intensities computed with plane-parallel and spherically-symmetric Atlas models (characterized by the three fundamental 
parameters L t , M+ and R+) are used to reexamine the existence of the fixed yu-point for the parametrized curves. 
Results. We find that the intensities from our spherical models do not have a fixed point, although the curves do have a minimum 
spread at a yii-value similar to the parametrized curves. We also find that the parametrized curves have two fixed points, /j t and /j 2 , 
although jj, 2 is so close to the edge of the disk that it is missed using plane-parallel atmospheres. We also find that the spherically- 
symmetric models appear to agree better with published microlensing observations relative to plane-parallel models. 
Conclusions. The intensity fixed point results from the choice of the parametrization used to represent the limb darkening and from 
the correlation of the coefficients of the parametrization, which is a consequence of their dependence on the angular moments of 
the intensity. For spherical atmospheres, the coefficients depend on the three fundamental parameters of the atmospheres, meaning 
that limb-darkening laws contain information about stellar atmospheres. This suggests that limb-darkening parameterizations fit with 
spherically-symmetric model atmospheres are powerful tools for comparing to observations of red giant stars. 

Key words, stars: atmospheres / stars: late type 



1. Introduction 

The variation of the specific intensity over a star's disk, com- 
monly called limb darkening or center-to-limb variation, is an 
important function of the physical structure of a stellar atmo- 
sphere. Because it is difficult to observe limb darkening for 
stars other than the Sun, it is common to represent limb dark- 
ening by analytic expressions whose coefficients are determined 
by matching the I,\ip) from model stellar atmospheres. Here 
fi = cos 6, where 6 is the angle between the vertical direction at 
that point on the stellar disk and the direction toward the distant 
observer. The most basic form of the law is the linear version, 
depending on fi to the first power, such as 



hU* = 1) 



or, more generally, 



h(M = 1) 



1 -u(l -fi). 



(1) 



(2) 



Limb-darkening laws have subsequently become more complex 
by including terms with yu, or its alternative r = sinfl, raised 
to higher integer powers as w ell as both yffi and log(yu) terms 
(IClaretll2000t lHowarthll2010l) . In addition to minimizing the fit 



to the model's hip), the laws have also introduced various con- 
straints, such as enforcing flux conservation at each observed 
wavelength. The intricacies of limb-darkening laws have in- 
creased as model stellar atmospheres have advanced. 

IHevrovskvl d2000l 120031 120071) . motivated by the potential of 
gravitational microlensing to provide particularly detailed mea- 
surements of stellar limb darkening, investigated the optimum 
method of extracting limb darkening from the data. One curi- 
ous feature that emerged from these studies is the existence of 
a fixed yu location on the stellar disk where the fits to the nor- 
malized model intensities have nearly the same value indepen- 
dent of wavelength or the model's or \ogg. This is true for 
both t he principal-component analysis developed bv IHevrovskvl 
(2003) and the more traditional linear limb-darkening law used 
by others. 

iHevrovsk?! (120031) used models computed with the Atlas 12 
code dKuruczlll996l) to construct his fitting procedure. These 
models use detailed opacity sampling to include many tens of 
millions of spectral lines in the radiative trans fer, but they s till as- 
sumed plane-parallel geometry, even though IHevrovskvl (120031) 
targeted red giants with T e s < 4000 K and log g < 1.0. These 
stellar parameters are exactly those where the assumption of 
plane-parallel geom e try sh ould break down. In the study by 
IClaret & Hauschildtl (120031) . which did use the spherical mod- 
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els of Hauschildt et al. ( 1999), the focus was on models with 
logg > 3.5, for which spherical extension is minimal. A broader 
study of limb darkening and the fixed yu-point using spherically 
extended model atmospheres is clearly needed. 



2. SAtlas model atmospheres 



iLester & Neilsonl (|2008) have developed the SAtlas code, a 
spherically extended version of Atlas. This code shares with 
Atlas the properties of static pressure structure, LTE popula- 
tions and massive line blanketing represented by either opacity 
distribution functions or opacity sampling (the faster opacity dis- 
tribution function version is used here). The spherical program 
differs from the plane-parallel Atlas by allowing the gravity to 
vary with radial position, and by computing the radiative trans- 
fer along rays through th e atmosphere i n the direction of the dis- 
tant ob server using the iRvbickil (Il97ll) version of the iFeautrieJ 
(11964 method, which accounts for the radial variation of the an- 
gle between the vertical and the direction of the ray. The struc- 
ture of the atmosphere is computed using a total of 8 1 rays whose 
angular spacing is determined by two factors. The first is the 
distribution of rays chosen to represent the "core" of the atmo- 
sphere, the region where the lower boundary condition for ra- 
diative transfer is the diffusion approximation. [Tester & Neilsonl 
(2008) found that different distributions of the core rays had no 
effect on the structure of the atmosphere, and so elected to use 
equal spacings in yiz. The remainder of the rays are tangent to the 
atmospheric levels at the stellar radius perpendicular to the cen- 
tral ray toward the observer. These rays are projections of the 
radial spacing of the atmospheric levels, which is logarithmic in 
the Rosseland mean optical depth. Because this distribution of 
rays is set by calculating the structure of the atmosphere, it is 
not necessarily optimal for studying limb darkening. Therefore, 
after computing the structure of the atmosphere, the surface in- 
tensities are derived from the structure rays by cubic spline in- 
terpolation for any desired number of r ays w ith any desired dis- 
tribution over the disk. As iHevrovskvl (120071) has demonstrated, 
cubic spline provides an excellent interpolation method. We have 
created our surface intensities at 1000 points equally spaced over 
<n < 1. 



3. Fixed point of stellar limb-darkening laws 



iHevrovsk?! (120001) found the fixed point in the limb-darkening 
profiles of both the Sun and i n plane-parallel Atlas models of 
red giants. iFields et all (120031) found a similar fixed point in 
their analysis of a K3 giant using spherical Phoenix models, al- 
though they excluded the limb from their intensity analysis be- 
cause the low intensity near the limb contributed almost noth- 
ing to the observations they were analyzing. Their truncation 
point ranged from r = 0.998 for \ogg = 3.5, corresponding 
to = 0.063, to r — 0.88 for logg = 0.0, corresponding to 
fi = 0.475. However, the truncation eliminated that part of the 
surface brightn ess that deviates m ost strongly from the plane- 
parallel model. IFields etail J2003) stated that the fixed point is 
a generic feature of any single-parameter limb-darkening law 
that conserves flux. For example, if the limb darkening is rep- 
resented as f(/j.) = I(n)/F = 2[1 + Ax(jj)], the fixed point is 

/"Axed = / _1 [2 Jo f(p?)fj?d(i']. Although they also state that the 
fixed point is not required by the two-parameter law they em- 
ployed, they found a fixed point in the microlensing observations 
they were modeling using such a law. 



To explore the parametrization more closely, we begin with 
the sa me two-parameter n ormalized limb-darkening function 
used bv IFields etail d2003l) . 



where *H in the Eddington flux, defined as 



<H 



(3) 



(4) 



A fixed point requires, for any two arbitrary models, that the 
intensity profiles satisfy 



him) = him), 

or in terms of Eq. [3] 



1 - Ai 1 1 - -// 



(5) 



Bill- 



(6) 



For this to be true, A and B must be linearly dependent, A - 
aB + B or A\ — A2 = a(Bi - #2). Substituting this relation into 
Eq. |6]yields an equation for the fixed-point 



3 5 ^_ 

-Qr//o + - Y/Uo - 1 - a = 0. 



(7) 



We need to verify that A = f(B) and to understand the properties 
of the parameter a. 

Applying a general least-squares method to Eq. [3] 



X' 



-S 



(8) 



we determine the coefficients A and B from the constraint equa- 
tions 



dA £-1 



x|l--A/i|=0 



(9) 



and 



8B 



x 1 1 - 1 Vj^l) = 0. 



(10) 



Multiplying Eq. [9] and Eq. [10] by Afi and then converting the 
summation to integration, J d/i, we create the equations that 
determine A and B, 



J 1 1 
— - + -A + -B -1=0 
2'H 4 6 



and 

_ 5 
~ 4 



j 1 11 1 

+ — - + -A+ —B - - = 0, 

2'H 6 96 6 



(ID 



(12) 
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where J is the usual mean intensity. It is also useful to define an 
angular pseudo-moment of the intensity, 



P 



Jo 



/(/<) A/£d//, 



which allows Eq.Q~2]to be written as 



P 



2«H 



J 1 11 1 
+ — - + -A + —B - - = 0. 

IfH 6 96 6 



(13) 



(14) 



Both / and P are determined by the properties of the model stel- 
lar atmosphere or the star whose observations are being fit by a 
limb-darkening law. Equation Q~T] and Eq. [12] clearly show that 
the coefficients A and B are uniquely determined by J and P. 
However, for A and B to be linearly dependent, that is A = aB+B, 
J and P must also be linearly related. 

The relation of J to P can be understood by following th e 
discussion of the diffusion approximation in Mihalas] (Il978h . 
This begins by representing the source function by a power 
series, which leads to the intensity being given by 7 v (T v ,/f) = 
XT=oM" [d n B v (T)/dT" v ] (Mihalas 1978, Eq. 2-88). Using this ex- 
pansion in the definitions of J and P, and keeping just the first- 
order term, gives the familiar results that J v as B v and K v » B v /3 
plus the additional result that P v a 2B v /3. Eliminating B v be- 
tween the J y and P v expressions gives 



P v = 2J v /3, 



(15) 



confirming the desired correlation. Of course, as the atmosphere 
thins out toward space the diffusion approximation becomes less 
accurate, but, as will be seen, the basic correlation of P and J is 
still there. 

Using Eq.[l5]we can combine Eq.fTTIand Eq.Q~2]to find that 



-0.694B. 



(16) 



In terms of the notation used earlier, a — -0.694 and B — 0, or 
equivalently AA = -0.694AB. Using this value for a in Eq. [7] 
leads to the equation for the fixed point, 



1.084/45 - 0.925/zq + 0.094 = 0. 



(17) 



This quadratic equation yields two fixed points, not one. The 

o 

solutions are yL\ = 0.736, corresponding to 0\ = 42.61 and 

o 

r\ - 0.677, and pL 2 = 0.118, corresponding to 6 2 = 83.22 and 
r 2 = 0.993. 

The results above used the diffusion approximation to es- 
tablish the relationship between J and P. To generalize these 
results, we define a new variable, 77, as 



= / 7V/7d/f _ p 

11 = 7 ~ 7' 



(18) 



Note that 77 is a generalization of the relation in Eq.[l5] Using the 
definition of 77, we replace the variable P in Eq. [14] by rjj, and 
then Eq.Q~T]is use d to eliminate the J jlfH terms. The resulting 
equation is rearranged to the form A = /(B), which is a linear 
equation of the same form as the previous equation for A, namely 
A = aB + B. Comparing the two equations for A we identify 



a = —- 



(5/7/4- l)/6+ 11/96 
(5/7/4 -l)/4+ 1/6 ' 



(19) 



Equation[l9]enables us to explore how a changes with rj. In ad- 
dition, because EqQis a function of a only, we are also able to 
determine the dependence of the two fixed points, //[ and jU2, on 
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Fig. 1. Top panel: dependence of the slope of the function A = 
aB +Bon the variable ij = J~ 1 fl ^dfi. Bottom panel: variation 
of the fixed points with r\. 



r\. These dependencies are shown in Fig. [T] Because the value 
of 77 can be set to values other than 2/3, the assumption of the 
diffusion approximation is no longer being used. 

Figure [T] shows that, except for rj m 0.267, the value of a 
is in the narrow range -0.72 < a < -0.68, which leads to the 
fixed points being 0.72 < m < 0.74 and 0.08 < fi 2 < 0.14. The 
divergence of a to +00 as 77 approaches as 0.267 is due to the 
denominator of Eq.[T9lapproaching zero as 77 — > 4/15. However, 
77 can never equal 4/15. This is a consequence of the following 
inequalities, which are true because fi < 1 : 



J = j Idfi > j I ^d/j > j Ij/dp = K. 

Dividing through by J this becomes 
K 

1 > 77 > — . 



(20) 



(21) 



Deep in the atmosphere, where the diffusion approximation 
holds, K/J = 1/3. Moving toward the surface the atmo- 
sphere becomes more transparent and the radiation becomes less 
isotropic, with the consequence that K/J > 1/3. As a result, 



77 > 1/3 > 4/15 = 0.267, 



(22) 



and the divergence cannot occur. Equation[T9]can be used to set 
bounds on a by using the lower bound on 77 from Eq. [22]and the 
upper bound on 77 from Eq.l2TI Theresultis -5/6 < a < -15/22. 
Using this range of a in Eq. [7] we find that the ranges for the 
fixed points are 0.708 < fii < 0.741 and 0.025 < /j 2 < 0.130, 
which are consistent the values for //; and fi 2 found previously. 
This indicates that the fixed points exist in the limb-darkening 
relations and are stable to differences in the intensity profile of 
the stellar atmosphere. It also suggests that the fixed points are 
not a result of a particular dominant opacity source. The only 
requirement is that 77 > 1/3. 
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Fig. 2. Wavelength dependence of the mean values of (a) the ra- 
tio of the pseudo-moment to the mean intensity, 77, (b) the value 
of a that defines the fixed points, (c) the first and (d) second fixed 
point for the grid of plane-parallel Atlas model atmospheres. 
The error bars represent 1-cr deviation from the mean. 



We verify this analytic result by computing a grid of Atlas 
plane-parallel model stellar atmospheres with 3000 K < T e g < 
8000 K and -2 < \ogg < 3. For this grid of models, the 
mean value of 77 in the V-band is fj pp = 0.716 ± 0.004 which 
leads to a pp = 0.691 + 0.001, m = 0.7375 ± 0.0002 and 
l-i2 = 0.1200 + 0.0003. This is very similar to what is found as- 
suming the diffusion approximation. The difference between the 
model V-band 77 and the value of 77 from the diffusion approx- 
imation is due to the fact that plane-parallel models do not en- 
force the diffusion approximation at all depths in the atmosphere. 
To determine the dependence on wavelength, we used the plane- 
parallel grid of atmospheres to determine the mean value and 1- 
cr standard deviation of 77, a and the two fixed points fii and /i2, 
with the results shown in Fig. [2] The values of the fixed points 
clearly do not vary significantly as a function of wavelength nor 
as a function of effective temperature and gravity. Furthermore, 
it must be noted that the value of 7/ approaches 2/3 as A — > 00. 

These results can be generalized to show that the fixed point 
occurs in other limb-darkening laws, such as a quadratic law 
where 7/2'TY = 1 - A(\ — jj.) — B(\ - 2fi 2 ). Repeating the deriva- 
tion above, the pseudo-moment, V, is replaced by K, a higher 
angular moment of the intensity. Because the diffusion approx- 
imation analysis also yields a linear relation between J and K, 
J = 3K, we again find that A = f(B). For a more general law, 
such as I/2<H = 1 - A(l - fx) - B[\ - \(n + 2)/j"], we would 
find a similar connection between the coefficients A and B. We 
would even find fixed points for any law that is a combination 
of a linear term and any function that can be represented by a 
power-la w series, such as e^, s i milar to the limb-darkening law 
tested by IClaret & Hauschildl (120031) . Therefore, we conclude 
that the fixed points are a general property of this family of limb- 
darkening laws. 



Fig. 3. Models computed wit h the SAtlas code p lotted over the 
stellar evolution tracks from Girardi et al. l (5000l) . 

4. Limb darkening and SAtlas model atmospheres 

To conduct a broader survey of limb darkening, we have com- 
puted several larger cubes of solar-composition model atmo- 
spheres, increasing the range of luminosity, mass and radius as 
well as using micro turbulent velocities of 0, 2 and 4 km/s. The 
luminosities, masses and radii of these models cover a significant 
portion of the Hertzsprung-Russell diagram where stars have ex- 
tended atmospheres, as shown in Fig. [3] which also shows the 
evolutionary tracks of Girard fet al.l(l2000l) for comparison. The 
number of models in each cube is large. For example, there were 
2101 models in the cube with v tur b = 2 km/s and the cubes for 
the other values of v tur b have similar numbers of models. For 
each model in each cube we computed the intensity at 1000 
equally spaced /u points at all wavelengths from the far ultra- 
violet to the far infrared; the results are shown in the left column 
of Fig. |4] for the B,V,R,I and //-bands. These intensity curves 
have similar be havior as the Phoenix mod el atmosphere inten- 
sity curves that lOrosz & Hauschildtl (2000) studied. We also fit 
each of the model limb-darkening curves with the linear-plus- 
square-root limb darkening law given by Eq.|3] and this is shown 
in the right column of Fig. |4] Note that the A and B coefficients 
of this parametrization have been shown in Eq. QT|and Eq. fT2lto 
be functions of J and P, not the structure of the intensity profile. 
The results for the v tul -b = and 4 km/s cubes of models are very 
similar. This is not surprising; in the Atlas codes, the turbulent 
pressure is P tu ,b = pv\ mh , indicating that the turbulent velocity 
has a similar effect on the temperature structure of a stellar at- 
mosphere as the gravity, that is a change of turbulent velocity is 
equivalent to a change o f gravity, similar to the discussion from 
iGustafsson et all (120081) for the MARCS code. In terms of the 
best-fit relations, a change of turbulent velocity causes a small 
change for the coefficients in the same manner as a change in 
gravity for the same effective temperature and mass. Therefore, 
we can just explore the grid of models atmospheres for only one 
value of v tur b. The cubes of models used in this survey cover a 
much wider range of atmospheric parameters than was used in 
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fi IL 

Fig. 4. Left panel: Limb darkening, characterized by the normalized surface intensity, 7(yu)/2'H, computed for the cube of 2102 
model atmospheres for v tul -b = 2 km/s for the BVRI and //-bands. Right panel: Fit to these same intensities using the limb-darkening 
law I(jj)/2'J-( = 1 — A(l - - B(\ — | ^fjT) for the same wavebands. 



the only pr evious investigation using s pherical model stellar at- 
mospheres (IClaret & Hauschi ldt 2003), which also truncated the 
profiles at the limb to achieve better fits. 

The model intensities in the left column of Fig.|4]are always 
positive, but the intensities in some of the parametrized fits in 
the right column become slightly negative toward the limb. This 
does not happen for plane-parallel model atmospheres because 
the slope of the intensity profile {dl/d/j. or dl/dr) has a constant 
sign. For example, using /i as the independent variable the slope 
of the intensity profile is always positive, varying from zero at 
the center of the disk to oo as /j — > 0. For spherical atmospheres, 
the center-to-limb variation changes sign because of a slight in- 
flection, which is apparent in the curves of the left column of 
Fig. [4] The more complex intensity profile and the use of an 
equal-weighting^ 2 -fit (i.e. the same number of /v-points near the 
center as near the limb) leads to some slightly negative intensi- 
ties in some best-fit relations. 



Figure |4] shows that the actual limb darkening com- 
puted from the spherical models is more complex than the 
parametrized representation, and it is obvious that the analysis of 
specific observations that contain intensity information should 
be cautious about using the parametrized fits. Although fits to 
plane-parallel models appear to be very good, the models are 
less realistic physically than the spherical models. Therefore, 
using what seems like a well-fitting law may introduce hidden 
errors that could compromise the conclusions of the analysis. 
The quantities derived from observations of binary eclipses and 
planetary transits may be more uncerta i n than had been thought 
(see, for example , iKnutson et al.ll2.Q07t ISouthworth et al.l [2007: 
ISouthworthl2"008h . 

Using the same grid of spherical Atlas models used in Fig. [4] 
we computed the mean values and 1-cr standard deviations of rj, 
a and the fixed points as a function of wavelength. The results 
are shown in Fig. [5] which can be compared with the plane- 
parallel results shown in Fig. [2] There are two key differences 
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Fig. 5. Wavelength dependence of the mean values of (a) the ra- 
tio of the pseudo-moment to the mean intensity, 77, (b) the value 
of a that defines the fixed points, (c) the first and (d) second fixed 
point for the grid of spherically symmetric Atlas model atmo- 
spheres. The error bars represent l-cr deviation from the mean. 



between the plane-parallel and spherical models. The first is that 
the mean values of 77 are larger for the spherical models, which 
leads to the differences between the mean values of a and the two 
fixed points. The second difference is that the standard deviation 
for 77 as a function of wavelength is much larger for the spherical 
models than for the plane-parallel models. This suggests that 77 
varies much more as a function of effective temperature, gravity 
and mass. Clearly, the fixed point is much more constrained in 
plane-parallel model stellar atmospheres. 

A specific example of an appli cation of limb dark ening to 
an observation is the analysis by iFields etall (120031) of mi- 
crolensing observations. In Fig. [6] we compare the V-band limb- 
darkening relations from our cube of spherical models with 
the V -band limb-darkening relation determined by Fields et al. 
(2003) from observations. It is clear that the Fields et al. (2003) 
V-band limb-darkening law agrees well with the limb-darkening 
from the spherical models, but, not as well with the plane- 
parallel models. Furthermore, although the curves for the model 
stellar atmospheres narrow to a waste rather than to a point, 
the location of minimum spread coincides with the location of 
the observed /i-position in the V-band limb-darkening relation 
within the uncertainty of the observed coefficients. Even though 
the uncertainty of limb-darkening relations from microlensing 
observations is large, these results suggest that the IFields et aU 
(2003) limb-darkening observations are probing the extended the 
atmosphere of the lensed red giant star. We conclude that the 
parametrized laws, although more simplified than the computed 
limb-darkening curves, are useful for understanding how limb 
darkening depends on the fundamental properties of the stellar 
atmosphere. 




1 0.8 0.6 0.4 0,2 



Fig. 6. Co mparison of t he V-b and limb-darkening relation deter- 
mined by Fiel ds et all ((2003) from observations (black dashed 
curves with errorbars) with the limb-darkening relations com- 
puted from SAtlas model stellar atmospheres (red solid curves) 
and relations computed from Atlas model stellar atmospheres 
(blue dotted curves), represented using the same linear-plus- 
square-root law. The insert shows a magnified view of the region 
of the fixed /./-point. 



To explore further the information content of the limb dark- 
ening, we recall that Eq. [16] showed that the coefficients A and 
B of the parametrization are linearly correlated. In Fig. [7] we 
plot the limb-darkening coefficients A and B for the V, /, and 
//-bands computed with both the plane-parallel Atlas and the 
spherical SAtlas models. For each wavelength th e erro r box 
shows the range of A and B found by Fiel ds et alj (|2003) from 
their observations. For the V-band, the spherical limb-darkening 
coefficients overlap with the observed fit, but the plane-parallel 
coefficients do not, while for //-band both plane-parallel and 
spherical models agree with observations. The microlensing 
observations at the longer wavelength might not be sensitive 
enough to probe the low intensity limb of the star, making the 
star appear consistent with limb-darkening of the plane-parallel 
model atmospheres. In the /-band neither set of models agree 
with the o bservations, but the /-band data provide a weak con- 
straint. As IFields et al.l (12003) noted, the /-band time-series mi- 
crolensing observations were a composite from multiple sites, 
and removing data from any one site changed the results signifi- 
cantly. 

In the V-band, the limb-darkening coefficients from 15 
model atmospheres fall within the observational box. These 
models have log g = 2.25 - 3 and r eff = 3400 - 3600 K. 
For the //-band data there are four models within the observa- 
tional box; these have same range of gravities but are slightly 
cooler, r e ff = 3000 - 3100 K. It is interesting that the models 
that agree with the observatio ns are those with gr avitie s that are 
consis tent with the results of IFields et al.l d2003l) and lAn et al.l 
d2002l) . This suggests that observations using the limb-darkening 
parametrization can probe the spherical extension of stellar at- 
mosphere via the fixed point, fi\, and potentially probe the fun- 
damental parameters of stars via the limb-darkening coefficients. 

Using a cube of models with a given value of v tul -b, we exam- 
ine the dependence of the primary fixed point, yUi , on the effective 
temperature and surface gravity. Note again that there are three 
basic parameters characterizing the spherical atmospheres: L*, 
M+ and /?*. This means that values of T e g and logg are degen- 
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Fig. 8. Left: As a function of surface gravity, the top panel shows the value of the primary fixed point, [i\ , for the linear plus square 
root parametrization, for spherical atmospheres of varying effective temperature. The bottom panel shows the dependence of the 
normalized V-band intensity of the fixed point. Red hilled circles are r e fr = 3000 K, green open squares 4000 K, blue open circles 
5000 K, magenta downward pointing triangles 6000 K and pale blue upward pointing triangles are 7000 K. The black crosses show 
the behavior of the fixed point from the grid of plane-parallel model atmospheres for comparison Right: As a function of effective 
temperature, the top panel shows the fixed point, [i\ , and the bottom panel shows the dependence of the V-band intensity of the fixed 
point. Red triangles represents models with log g — 0, green squares log g = 1, and blue circles are log g — 2. Again, black crosses 
represent the fixed point from plane-parallel model atmospheres. 



erate, but they are easier to show on a two-dimensional surface. 
In the top left panel of Fig. [8] we plot the value of the primary 
fixed point, fi i , as a function of log g. At each value of log g there 
are values of T e ff ranging from 3000 to 7000 K, although there 
can be multiple values because of the parameter degeneracy. In 
the bottom left panel of Fig.[8]we plot Iy(/ui)/2'J-ly, the normal- 
ized intensity of the fixed point in the V-band. For each surface 
gravity the values of /yi and the normalized intensity at the fixed 
point both show a steady progression as T e s changes, except for 
r eff = 3000 K. We suspect that the behavior for T eB = 3000 K 
is due to a change in the dominant opacity source for our coolest 
models, possibly water vapor. Models with T e ff > 3000 K have 
H~ as the dominant continuous opacity, but in the coolest mod- 
els there are fewer free electrons available to form H , and the 
formation of H2 reduces the pool of hydrogen atoms. We also 
plot the fixed point and intensity at the fixed point from plane- 
parallel model atmospheres for comparison. It is clear that the 
fixed point from fits to plane-parallel models varies much less 
than fits to spherical models. 



On the right side of Fig. [8] we reverse the parameters and plot 
the dependence of the primary fixed point and the normalized in- 
tensity of the fixed point as a function of T e g. At each value of 
the effective temperature, the values of log g are 0, 1 and 2 in cgs 
units. In the top right panel of Fig. [8]we see that there is essen- 
tially no variation of the value of [i\ with T e s for all three surface 
gravities until the lowest effective temperature is reached. There 
is an obvious displacement of the value of /ii for each gravity and 
also a spread in /ii because of parameter degeneracy. However, 
for T e ff < 3500 K the value of /u\ drops for all surface gravities. 
The bottom right panel shows that the normalized V-band inten- 
sity of the fixed point shows a similar behavior. For T e ff > 3500 
K there is little variation with r e ff, but there is an offset and a 
spread that depends on the surface gravity. Cooler than 3500 K 
the value of the normalized intensity drops noticeably. 

It is clear that the fixed point [i\ and the intensity at the fixed 
point l(yi\)j2'H are functions of gravity and effective temper- 
ature for spherically-symmetric models and are roughly con- 
stant for plane-parallel model atmospheres. This indicates that 
the best-fit coefficients of the limb-darkening law vary mostly 
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Fig. 7. Red plus symbols represent the limb-darkening coeffi- 
cients A and B computed from SAtlas model intensity profiles 
for the V, I, and //-bands, and green open squares represent 
the coefficients computed with plane-parallel Atlas models. The 
box in each plot shows the ra nge of the coefficien ts derived from 
microlensing observations by Fiel ds et al.l(l2003l) . 

because of the geometry of the model atmosphere. A spherical 
model atmosphere predicts a smaller intensity near the limb of 
the stellar disk relative to a plane-parallel model with the same 
effective temperature and gravity. To predict the same emergent 
flux, the intensity must be larger at the center of the disk, hence 
the temperature at the base of the atmosphere must also be larger 
for the spherical model. Therefore, the temperature structure of 
the atmosphere also varies. However, this is a secondary effect 
and the geometry of the atmosphere is most important in de- 
termining the value of the /yi and /( J t/i)/2'K. The geometry of 
spherical models leads to smaller values of the pseudo-moment 
and the mean intensity because the intensity is more centrally 
concentrated. This suggests that n\ and /(yUi )/2"7/ depend on the 
atmospheric extension. We will explore how the fixed point and 
intensity relate to the extension and fundamental stellar proper- 
ties in greater detail in a future article. 

5. Conclusions 

We have explored limb darkening using large cubes of spheri- 
cal stellar atmospheres spanning the parameters L*, M* and /?* 
covering the cool, luminous quadrant of the Hertzsprung-Russell 
diagrams (Fig. [3). These models have also used three different 



values of the microturbulent velocity. For each model, the center- 
to-limb variation of the surface intensity has been calculated at 
1000 equally spaced fi values spanning the range from 1 to for 
every wavelength used to compute the model structure. 

Parametrizing the center-to-limb variation with a flux- 
conserving linear-pl us-square-root limb -dar kening law, we con- 
firm the findings of iHevrovskvl (l2000l) and iFields et al.1 (l2003h 
that there is a fixed [i\ point through which all the intensity 
curves pass. However, when we plot the surface intensities di- 
rectly, without using a fitting law, there is no fixed point, al- 
though the distribution of curves does narrow to a waist close 
to the same value of yUi (Fig.|4]i. 

The apparent fixed point is a result of the least-squares fit- 
ting procedure where the two parameters of the law are depen- 
dent on two properties of the stellar atmosphere, the mean inten- 
sity, J, and the pseudo-moment, J I ^dfi. For the temperature 
range 4000 - 8000 K, the mean intensity is correlated with the 
pseudo-moment, which means that the two coefficients are also 
correlated, leading to the existence of the fixed point. 

The lack of a well-defined fixed point in the surface intensity 
distribution for spherical model atmospheres suggests that the 
three fundamental parameters of the atmospheres affect the limb 
darkening in a way that is not encountered in the two-parameter 
plane parallel atmospheres. 
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